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Wc present a numerical study of lasing modes in diffusive random media with local pumping. The reab- 
sorption of emitted light suppresses the feedback from the unpumped part of the sample and effectively 
reduces the system size. The lasing modes are dramatically different from the quasimodes of the passive 
system (without gain or absorption). Even if all the quasimodes of a passive diffusive system are extended 
across the entire sample, the lasing modes are still confined in the pumped volume with an exponential 
tail outside it. The reduction of effective system volume by absorption broadens the distribution of decay 
rates of quasimodes and facilitates the occurrence of discrete lasing peaks. © 2008 Optical Society of America 
OCIS codes: 290.1990,140.3460,290.4210 



Over the past few years, there have been many studies 
on random lasers with coherent feedback. 1 In a strongly 
scattering active medium, recurrent scattering events 
could provide resonant feedback for lasing. With suffi- 
cient gain, lasing oscillation might occur at discrete fre- 
quencies that are determined by the interference of scat- 
tered light. In our initial experiments with highly disor- 
dered semiconductor powder and polycrystallinc films, 
the transport mean free path was short, leading to small 
laser cavities. 2 Our experimental observation was well re- 
produced in the numerical simulation of lasing in local- 
ized modes. 3-5 However, lasing with coherent feedback 
was realized also in weakly scattering random media. 6-8 
Tight focusing of pump light was necessary to observe 
discrete lasing peaks, namely, the pump beam was fo- 
cused to a region of size much smaller than the entire 
sample. Imaging of laser light on the sample surface re- 
vealed that the lasing modes were not extended over the 
entire random medium, instead they were located in- 
side the pumped region with an exponential tail outside 
it. 9 Since the quasimodes of a random system far from 
the onset of localization were usually extended states, 
the lasing modes were regarded as some types of anoma- 
lously localized states, either almost localized states, 10, 11 
or prelocalized states. 9 ' 12 

The anomalously localized states should be rare in the 
diffusive samples. Yet no matter where on the sample the 
pump beam is focused, we always observe lasing modes 
that are spatially confined in the pumped region. More- 
over, the lasing threshold does not fluctuate much as 
the pump spot is moved across the random medium. 
We believe the contradiction to the theory of anoma- 
lously localized states originates from the assumption 
that lasing occurs in the quasimodes of the passive ran- 
dom medium. In this letter, we will demonstrate with 
numerical simulation that this assumption is not valid 



when absorption at the emission wavelength is significant 
outside the pumped volume. The reabsorption of emitted 
light suppresses the feedback from the unpumped part of 
the sample and effectively reduces the system size. The 
lasing modes are formed by the interplay between gain 
and absorption, and therefore, are dramatically different 
from the quasimodes of the passive system (without gain 
or absorption). Even if all the quasimodes of a passive 
diffusive system are extended across the entire sample, 
the lasing modes arc still confined in the pumped volume 
with an exponential tail outside it. 

We use the finite-difference time-domain (FDTD) 
method to simulate lasing in the transverse magnetic 
(TM) modes of two-dimensional (2D) random media. 
The disordered system is a collection of dielectric cylin- 
ders placed at random in vacuum. The diameter of the 
cylinders is 160nm, the refractive index n = 2, and the 
filling fraction is about 50%. The total size of the sys- 
tem is 9.2/xm x 9.2/im. The random medium is sur- 
rounded by vacuum, which is terminated by uniaxial 
perfectly matched absorbing layer. 13 The wavelength of 
interest is around 650nm. To verify light transport in 
the 2D system is diffusive, we calculate the continuous 
wave (CW) transmission through a slab of the random 
medium with the thickness L = 9.2/im. The ensemble- 
averaged intensity profile (I(z)) is fitted well with the 
diffusion expression (L + z$ — z)/(L + 2zo), 14 where z 
is the coordinate normal to the slab. The extrapolation 
length zq is related to the transport mean free path £ 
as Zo/i = 7r/4(l + R)/(l — R). R is angularly-averaged 
boundary reflection coefficient, which accounts for in- 
terface effects. 14 Using the expression for R in the 2D 
system, 15 we estimate z /£ ~ 1.85. The value of zq is ex- 
tracted from the fitting of (I(z)), that gives I ~ 1.3/j.m. 
The small value of IjL ~ 0.12 as well as the large value 
of k£ ~ 13 ensures the 2D random system of lateral di- 




Fig. 1. Spatial intensity distribution of (a) the quasi- 
modc with the longest lifetime in a passive diffusive sys- 
tem; (b) the (first) lasing mode with gain inside the cir- 
cular region near the center and no absorption outside it; 
(c) the (first) lasing mode with gain inside the circular 
region near the center and absorption outside it; (d) the 
(first) lasing mode with random medium beyond one l v 
(dashed circle) removed. 
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Fig. 2. Radial dependence of the angularly-integrated 
intensities of the modes (a,c,d) shown in Fig. ^ ( a ) is 
compared to the diffusive mode profile; the tail of (c) is 
fitted with an exponential decay; the falloff of (d) outside 
the random medium is compared to r _1 . 

mension L — 9.2/^m is in the diffusive regime. 

We start with a study of the quasimodes (eigen- 
modes of the Maxwell's equations) in the passive ran- 
dom medium. A short excitation pulse, whose spectrum 
is centered at Ao — 650nm, is launched in the center 
of the sample. After the excitation pulse is gone, the 
total electromagnetic energy stored inside the random 
system U (t) exhibits a non-exponential decay in time as 
a result of multi-mode excitation. However, after a suffi- 



ciently long time, U(t) changes to an exponential decay 
because only one mode with the longest lifetime is left 
inside the system. The exponential decay rate of U (t) is 
equal to the decay rate 7 of this quasimode. For example, 
for one realization shown in Fig.^a), the Fourier trans- 
form of electric field E(t) gives the wavelength of the 
quasimode A = 646nm. The spatial profile of the electric 
field becomes stable in time, and it reflects the wave- 
function of this longest-lived quasimode. As shown in 
Fig.^a), this quasimode is extended across the sample. 
When we remove some cylinders at the sample bound- 
ary, e.g., the cylinders outside a circle of radius L/2 from 
the mode center, the mode profile changes dramatically. 
The sensitivity of the quasimode to the boundary con- 
firms that it is an extended state. Fig.|21shows the radial 
dependence of the angularly- integrated intensity, I r (r) — 
j \E(r, 9)\ 2 rd8, where r is the radial coordinate. 6 is the 
polar angle, E(r, 9) is the electric field distribution of the 
quasimode. For comparison, the radial profile of the low- 
est diffusion mode cos [irx/ (L + 2zo)] cos [iry/ (L + 2zq)] 
is also plotted in Fig. [3 It describes I r (r) relatively well. 

Now we introduce gain and reabsorption into the 
random medium. In the previous studies, 3,4 the gain 
medium is modeled as four-level atomic system where 
the lasing transition is from the third level to the sec- 
ond level. In the absence of pumping, all the atoms are 
assumed to be in the first level. Since the electronic pop- 
ulation in the second level is zero, there is no absorp- 
tion at the lasing wavelength in the unpumped region. 
To simulate spatially non-uniform gain and reabsorption 
of the laser emission, we use the semi-classical Lorentz 
model. 13 The linear gain/absorption is modeled by nega- 
tive/positive conductance. 2 By introducing negative con- 
ductance to the pumped region and positive conductance 
to the unpumped region, we are able to describe both 
light amplification inside the pumped region and reab- 
sorption of the emitted light outside the pumped region. 
More specifically, the cylinders have the conductance 13 



1 + i(u) — ujo)T 2 1 + i(uj + ojq)T 2 

The sign of ao determines light is amplified or ab- 
sorbed, whereas the amplitude of ao sets the magnitude 
of gain/absorption, luq and I/T2 determine the center 
frequency and width of the gain/absorption spectrum, 
respectively. The absence of gain saturation in Eq. ^ is 
not crucial in our simulation, as our goal is to find the 
first lasing mode at or slightly above the lasing thresh- 
old. A seed pulse is launched at t = to initiate the 
amplification process. The lasing threshold is defined by 
the minimum gain coefficient (— ao) at which the electro- 
magnetic energy stored inside the random system grows 
in time and eventually diverges. 

We first consider the case of local pumping with no 
absorption outside the pumped region. We use the same 
random sample as in Fig.^a), and introduce gain to the 
central part marked by the circle in Fig. 1(b). The gain 
spectrum is centered at 650nm and has a width of 52nm. 



2 



The lasing mode shown in Fig. IHb) is identical to the 
quasimode of the passive system in Fig.^a). Although 
optical gain is concentrated within the circle near the 
center, the lasing mode is extended throughout the en- 
tire sample. As we reduce the pump area by decreasing 
the radius of the circle, the lasing threshold is increased, 
but the lasing mode profile remains the same. This re- 
sult indicates that the lasing mode in a diffusive random 
medium is the extended quasimode of the passive sys- 
tem, even when the pumped region is smaller than the 
mode size. 

However, the above statement is valid only when there 
is no absorption outside the pumped region, which is 
not the case in most experiments. For example, Rho- 
damine dye, which is widely used to provide gain for 
random lasers, has significant overlap between its ab- 
sorption band and emission band. Therefore, photons 
that are emitted by the excited rhodamine molecules in- 
side the pumped region may diffuse into the surround- 
ing unpumped region and be absorbed by the rhodamine 
molecules there. The absorption reduces the probability 
of light returning to the pumped region, thus suppresses 
the feedback from the unpumped region. To simulate 
the reabsorption, we introduce absorption outside the 
pump area. The bulk absorption length (without scat- 
tering) l a = 1.06/im. Fig. (He) shows the lasing mode 
profile, which is very different from that in Fig. [Ja) 
or^b). The wavelength of the lasing mode also differs 
by about 4nm. Therefore, the lasing mode in the pres- 
ence of reabsorption is a new mode, completely different 
from the quasimode of the passive system. Due to re- 
absorption outside the pump area, the lasing mode is 
confined more or less inside the pumped region. The ra- 
dial profile of the lasing mode, shown in Fig. [5J features 
an exponential decay outside the pumped region. The 
decay length, ~ 0.8/im, is equal to the diffusive absorp- 
tion length £ p = {lal/2) 1 / 2 . To confirm this result, we 
simulate transmission of CW plane wave through a slab 
of disordered absorbing medium, and obtain the same 
attenuation length. 

Performing calculations similar to those shown in Figs. 
IH2l for many disorder configurations, we find that the ob- 
served effect is indeed general. The reabsorption, which 
suppresses the feedback from the unpumped part of the 
sample, effectively reduces the system size. To check this 
conjecture, we remove all the random medium beyond 
one t p from the pump area [dashed circle in Fig. Q^d)], 
and repeat the calculation. The frequency and spatial 
profile of the lasing mode remain the same [Fig.QJd)], de- 
spite the drastic change of the random system. The radial 
profile of the lasing mode outside the random medium 
(Fig. [2J exhibits a trivial dependence. This result 
indicates the lasing mode is an extended state within 
the effective volume, V e f / , that is made of the optically 
active region plus the buffer layer of thickness l v around 
it. 

The reduction of the effective system volume leads to 
a decrease of the Thouless number S = 8v/Av, where 



5v and Av are the average mode linewidth and spac- 
ing respectively. In a 3D diffusive system bv oc V e ^ 

and Av cx therefore, S tx vHj. The smaller 

the value of 5, the larger the fluctuation of the decay 
rates 7 of the quasimodes. 12, 16 The variance of the de- 
cay rates 12 er 2 = (j) 2 /6, where the average decay rate 

(7} ~ D/V^j-f. We believe the broadening of the de- 
cay rate distribution along with the decrease of the total 
number of quasimodes (within V e ff) is responsible for 
the observation of discrete lasing peaks in the tight focus- 
ing condition. Despite its value is reduced, the effective 
Thouless number is still much larger than 1 due to weak 
scattering. As a result, the lasing modes are the extended 
states within the effective volume. Because cr 7 /(7) <C 1, 
the minimum decay rate is still close to (7}, leading to 
relatively small fluctuation of lasing threshold. 17-19 
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